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Cosmic microwave background (CMB) observations suggest the possibility of an extra dark radia- 
tion component, while the current evidence from big bang nucleosynthesis (BBN) is more ambiguous. 
Dark radiation from a decaying particle can affect these two processes differently. Early decays add 
an additional radiation component to both the CMB and BBN, while late decays can alter the 
radiation content seen in the CMB while having a negligible effect on BBN. Here we quantify this 
difference and explore the intermediate regime by examining particles decaying during BBN, i.e., 
particle lifetimes tx satisfying 0.1 sec < rx < 1000 sec. We calculate the change in the effective 
number of neutrino species, A^e//, as measured by the CMB, ANcmb, and the change in the effec- 
tive number of neutrino species as measured by BBN, ANbbn, as a function of the decaying particle 
initial energy density and lifetime, where ANbbn is defined in terms of the number of additional 
two-component neutrinos needed to produce the same change in the primordial *He abundance as 
our decaying particle. As expected, for short lifetimes (rx 0.1 sec), the particles decay before 
the onset of BBN, and ANcmb = ANbbn, while for long lifetimes {tx ^ 1000 sec), ANbbn is 
dominated by the energy density of the nonrelativistic particles before they decay, so that ANbbn 
remains nonzero and becomes independent of the particle lifetime. By varying both the particle en- 
ergy density and lifetime, one can obtain any desired combination of ANbbn and ANcmb, subject 
to the constraint that ANcmb > ANbbn- We present limits on the decaying particle parameters 
derived from observational constraints on ANcmb and ANbbn- 



I. INTRODUCTION 

Over the past decade, a "standard-model" cosmol- 
ogy has emerged, based, among other things, on preci- 
sion measurements of the fluctuations in the cosmic mi- 
crowave background (CMB) [l|, 0], observations of type 
la supernovae [1, and big bang nucleosynthesis (BBN) 
(seeRef. Q for a recent review). In the standard cosmo- 
logical model, the density of the universe is dominated 
at present by a cosmological constant (A) and cold dark 
matter (CDM), corresponding respectively to roughly 
70% and 25% of the total density, with the remaining 
5% in baryons. The present-day radiation content of the 
universe is negligible in comparison, but this radiation 
was the dominant component at early times. 

While the cosmological observations are generally con- 
sistent with this standard model, there are a few unre- 
solved issues. One of these involves the total radiation 
content of the universe. The energy density of the CMB is 
a simple function of the CMB temperature and is known 
to high accuracy. Similarly, given the observed number of 
light neutrinos {Ni, — 3), one can calculate the neutrino 
energy density. (In fact, the "effective" number of neutri- 
nos, Neff is sHghtly greater than 3 due to partial heating 
of the neutrinos in the early universe by electron-positron 
annihilation. Including these effects yields Ngff = 3.046 
S Q)- However, recent precision measurements of the 
CMB fluctuations are best fit by larger values of Neff- 
(For a discussion of the effect of N^f f on the CMB fluc- 
tuations, see Refs. [1, Q). The seven-year data from 
the Wilkinson Microwave Anisotropy Probe, combined 
with observations of baryon acoustic oscillations (BAO) 
and measurements of the Hubble parameter, Hq, give 
Neff = 4.34l°J^ (68% CL) [|. Observations by the At- 
acama Cosmology Telescope combined with BAO and Hq 



give Neff = 4.56 ± 0.75 (68% CL) Recent results 
from the South Pole Telescope combined with WMAP7, 
BAO, and Hq give N^ff = 3.86 ± 0.42 (68% CL) [HI- 
An analysis using comlained datasets in Ref. yields 
Neff = 4.08lJ^;^8 (95% CL). Since this additional radi- 
ation component cannot interact electromagnetically, it 
has been dubbed "dark radiation" . 

BBN is also quite sensitive to the total radiation con- 
tent in the universe, but here the evidence is more am- 
biguous. Recent calculations of the relic helium abun- 
dance by Izotov and Thuan [l^ and by Aver, Olive, and 
Skillman [T3| have reached opposite conclusions, with the 
former arguing for an additional dark radiation compo- 
nent and the latter concluding that the standard number 
of neutrinos suffices. A recent analysis by Mangano and 
Serpico [H gives ATV^// < 1 (95% CL). (See also the 
discussion in Ref. [l^). An overview combining CMB 
and BBN constraints on dark radiation can be found in 
Ref. while constraints on the physical properties of 
the dark radiation are discussed in Refs. [T^, [3 ■ 

Given these hints of new physics, a number of mod- 
els have been proposed to account for additional dark 
radiation. The simplest way to achieve this is to add ad- 
ditional relativistic relic particles, as suggested by, e.g., 
Refs. p^| - [2^ . In this case, the value for Neff determined 
by BBN and the CMB should be the same. Additional 
relativistic energy density can also be provided by a neu- 
trino chemical potential, as in the models discussed in 
Refs. mill. 

Alternately, if one wishes to produce an increase 
in Neff in the CMB, but retain the standard-model 
value for Neff in BBN, then an obvious possibility is 
the production of relativistic, non-electromagnetically- 
interacting particles from the decay of a massive relic 
particle after BBN. Such a scenario for the dark radia- 
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tion was considered by Ichikawa et al. [25|, and further 
elaborated by Fischler and Meyers [1^ and Hasenkamp 
[13] ■ If we let Nbbn denote the value of N^fj as mea- 
sured by BBN, and Ncmb be the value determined from 
the CMB and other low-redshift measurements, then 
these decaying particle models predict Nb bn ^ Ncm b , 
while additional stable relativistic degrees of freedom give 
Nbbn = Ncmb- 

In this paper, we fill in the gap between these two 
regimes by examining decaying particle scenarios in 
which the particle decays during BBN. Such models pro- 
duce a relation between Nbbn and Ncmb that varies 
from Nbbn = Ncmb at short particle lifetimes (when 
the particle decays before the onset of BBN) to Nb bn ^ 
Ncmb when the particle decays after the conclusion of 
BBN. In the next section, we give a detailed discussion 
of our calculation and present our results for Nbbn and 
Ncm s as a function of the decaying particle abundance 
and lifetime. Our conclusions, including observational 
limits, are discussed in Sec. III. 



II. DARK RADIATION FROM A DECAYING 
PARTICLE 

We assume a standard flat Friedman-Robertson- 
Walker model with the expansion rate given by: 
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where R is the scale factor and p is the energy density. 
To this standard cosmological model we add a nonrel- 
ativistic particle X, which is unstable and decays with 
lifetime tx- By assumption, X decays only into "invis- 
ible" relativistic decay products, which do not interact 
electromagnetically and can thus form the dark radia- 
tion (see Refs. [25l - l27l | for examples of such models). 
The effects of such decays during BBN were previously 
considered by the authors of Ref. [28|], and our treat- 
ment closely follows theirs. (For another early discussion 
of BBN with such decays, see Ref. [295.) The main focus 
of Ref. 2^ was the constraint from BBN that could be 
placed on these decaying particles, while in this paper, 
we will be interested in using such models to provide dark 
radiation, and to determine the relation between Nb b n 
and Ncmb as a function of the model parameters. 

We follow Ref. [S^ and parametrize the density of the 
decaying particle in terms of its number density relative 
to the entropy density, s, prior to decay {t tx)- 



where s is given by 
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In Eq. ([3]), T-y is the photon temperature, and is the 
effective number of "entropy" degrees of freedom, defined 



5«= E 9B{T,lT^f + {7/8) 9Fm/T^f, (4) 

bosons f ermions 

where the sum is over all relativistic bosons and fermions. 
In equation Q, gB and gp are the total number of bo- 
son and fermion spin degrees of freedom, and is the 
temperature of a given relativistic particle species. The 
important useful property of this parametrization is that 
Yx remains constant through the epoch of e+e^ annihila- 
tion, which occurs during the middle of BBN. During the 
BBN epoch, we have — 43/4, and nx/n^ — 19.36 Yx 
prior to e+e^ annihilation, while nx/n^ — 7.04 Yx after 
e^e~ annihilation. For comparison, the results of Ref. 
[2^ are parametrized in terms of nx jn^ prior to e^e~ 
annihilation. 

The equations governing the evolution of px and the 
decay-produced "invisible" radiation component, pdecj 
are [3l| 



dpx 
dt 

dpdec 

dt 

Equation ^ gives 



-iHpx - Px/tx, 
-AHpdec + Px/tx- 



(5) 
(6) 

(7) 



while equation ([B]) must be integrated numerically (al- 
though an analytic solution can be derived for t <ti tx 
[Slj). It is clear from Eqs. ([S])-® that pdec: and thus, 
the increase in N^ff, depends only on Yxmx and tx- 

Evolving these equations to calculate px and pdec, and 
converting pdec into an effective number of neutrinos, we 
derive ANcmb, the change in Ne/f as measured by the 
CMB, as a function of Yxmx and tx, where we confine 
our attention to the case where the X particle fully de- 
cays before last scattering. A contour plot of ANcmb as 
a function of Yxmx and tx is given in Fig. 1. 

In the limit where the decaying particles themselves 
never dominate the expansion, one can calculate ANc mb 
as a function of Yxmx and tx , as in Ref. |2^ . Rewriting 
the results of Ref. [28] in terms of our parameters, we 
find 



AN, 



CMB 



= 8.3 (Yxmx /MeY){Tx /secy Z"^. (8) 



A comparison of Eq. ([5]) with the results displayed in 
Fig. 1 shows that Eq. ([8]) is accurate to within ^ 10% 
for the curves displayed in Fig. 1. Note that our analytic 
expression for ANc mb differs significantly from that de- 
rived in Ref. [1^, as the latter used the "sudden decay 
approximation," in which all of the energy density of the 
decaying particle is taken to be converted into relativistic 
decay products at t = tx- 

Adding additional energy density during BBN affects 
all of the element abundances, but the effect is largest 
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FIG. 1: Contour plot of ANcmb, the change in the effective number of neutrinos determined by CMB observations due to a 
decaying particle with lifetime tx and energy density prior to decay parametrized in terms of Yxmx, where Yx is the initial 
number density of the particle relative to entropy density, and mx is the particle mass. Curves correspond to, from bottom to 
top, ANcMB = 0.1, 0.2, 0.3, 0.4, 0.7, 1.0, 1.5, 2.0. 



(relative to the accuracy with which the primordial abun- 
dances can be estimated) for ''He. Additional energy den- 
sity (either relativistic or nonrelativistic) increases the 
expansion rate, resulting in a larger relic neutron abun- 
dance, and increased "^He production [5|. Thus, the "^He 
abundance has long been used to constrain any addi- 
tional relativistic energy density, such as that produced 
by additional neutrinos [33 |. 

We use the Kawano version of the Wagoner [13, HI] 
Big Bang nucleosynthesis code to calculate the change in 
the primordial "^He abundance with the addition of px 
and pdec as given by Eqs. (O and (|6]), taking a baryon- 
photon ratio of 77 = 6.1 x 10"^^°. (Note that since we are 
presenting the change in the ^He abundance produced 
by the decaying particle, rather than the absolute he- 
lium abundance itself, our results are quite insensitive to 
the assumed value of 77; see, e.g.. Fig. 6 of Ref. Q). 
We determine the change in the number of relativistic 
neutrinos that gives exactly the same change in the ^He 
abundance as a decaying particle with a given abundance 
and lifetime. Thus, for any pair of values for Yxmx and 
Tx, we have a corresponding change ANbbn that pro- 
duces the same effect on BBN. (The change in the other 



element abundances can be ignored here). In Fig. 2, we 
give a contour plot of ANbbn as a function of Yxmx 
and Tx- 



III. CONCLUSIONS 

In comparing Figs. 1 and 2, we see that for tx ^ 0.1 
sec, ANbbn = ^Ncmb- In this short-lifetime limit, all 
of the decaying particle energy density is converted into 
dark radiation before BBN begins, so both BBN and the 
CMB "see" the same N^ff. 

In the opposite limit, tx ^ 1000 sec, the contours in 
Fig. 2 become horizontal lines. In this long-lifetime limit, 
all of the X particles decay after BBN, and the increase 
in the expansion rate that alters the *IIe abundance is 
due entirely to the energy density of the nonrelativistic 
particles before they decay. Thus, in this limit, ANbbn 
becomes a function only of Yxmx and is independent 
of Tx ■ Note that ANb bn never goes to zero in the long 
lifetime limit precisely because of this contribution to the 
expansion rate from the nonrelativistc particles. How- 
ever, by increasing rx, one can make ANbbn /^Ncmb 
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FIG. 2: Contour plot of ANbbn, the change in the effective number of neutrinos giving the same change in the primordial ''He 
abundance as a decaying particle with lifetime rx and energy density prior to decay parametrized in terms of Yxmx, where 
Yx is the initial number density of the particle relative to entropy density, and mx is the particle mass. Curves correspond to, 
from bottom to top, ANbbn = 0.1, 0.2, 0.3, 0.4, 0.7, 1.0, 1.5, 2.0. 



arbitrarily small. 

The interesting transitional regime, then, is precisely 
the one we have explored: 0.1 sec ^ tx ^ 1000 sec. This 
would be the regime of interest if more precise measure- 
ments of Nf,f f from the CMB and BBN yielded nonzero 
values for both ANcmb and ANbbn with ANbbn ^ 
ANcMB- It this case, it is possible to simply read off, 
from Figs. 1-2, values of Yxmx and tx that give the 
desired values for ANbbn and ANcmb- Note, however, 
that one always has ANcmb > ANbbn in this scenario, 
so this model can be falsified by observations contradict- 
ing this inequality. 

As an example, we show, in Fig. 3, the limits on Yxmx 
and Tx using the upper and lower bounds on ANcmb 
from Ref . [I4I and the upper bound on A b n from Ref. 

. As expected, the upper bound on ANb b n cuts into 
the region favored by ANcmb at short lifetimes, but cur- 
rent bounds are not sufficiently restrictive for this to be 
a major effect. Tighter bounds from future observational 
data will, of course, shrink this allowed region. 

Our results can be generalized to more complicated 



scenarios, such as a particle that decays into both dark 
radiation and electromagnetically-interacting particles. 
Even a small branching ratio into the latter can pro- 
duce a markedly different effect on BBN if the decay 
products are energetic enough to photofission the pri- 
mordial nuclei (see, e.g., Refs. j36l - l38l | and references 
therein). These effects, however, tend to be minimal 
for the lifetimes (rx ^ 10'^ sec) considered here, since 
the electromagnetically-interacting particles thermalize 
rapidly at high temperatures. Even in this case, how- 
ever, the electromagnetically-interacting decay products 
can heat the photons relative to the neutrino background, 
potentially decreasing N^f / instead of increasing it [3§| . 
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FIG. 3: Observational constraints from the CMB and BBN in the tx, Yxmx plane, where tx is the decaying particle lifetime, 
Yx is the initial number density of the particle relative to entropy density, and mx is the particle mass. The region between the 
two solid curves is allowed by the CMB bounds on A^e// from Ref. [l^ . while the area below the dashed curve is the allowed 
region from BBN limits on N^ff given by Ref. [l5| ]. The shaded (yellow) region satisfies both constraints. 
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